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ABSTRACT 

Brezin-Hikami contour-integral representation of exponential multidensities in finite A'^ Hermitian matrix model is a remarkable 
implication of the old Hermitian-Kontsevich duality. It is also a simplified version of Okounkov's formulas for the same multidensities in 
the cubic Kontsevich model and of Nekrasov calculus for LMNS integrals, a central piece of the modern studies of AGT relations. In this 
paper we use Brezin-Hikami representation to derive explicit expressions for the Harer-Zagier multidensities (from arXiv:0906.0036): 
^^J the only known exhaustive generating functions of all-genera Gaussian correlators which are fully calculable and expressed in terms 

of elementary functions. Using the Brezin-Hikami contour integrals, we rederive the 1-point function of Harcr and Zagier and the 
2-point arctangent function of arXiv:0906.0036. We also present (without a proof) the explicit expression for the 3-point function in 
terms of arctangents. Derivation of the 3-point and higher Harer-Zagier functions remains a challenging problem. 
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,3 ; 1 Introduction 

Gaussian correlators are among the central quantities in mathematical physics. For many decades they remain the 
subject of constant investigation and innumerable applications, still we are far from a satisfactory and exhaustive 
J^ ' description. The problem sounds very simple: to evaluate Gaussian integrals over N x N Hermitian matrices 

Cr''= ((tr (/.*)) =C, (2) 

C[°""= ((tr0Hr0^)) =a,-aCj (3) 

C^°r = ((tr 0Hr 0^"tr 0^)) = Q,k - Q,Ck - QkC, - C,kQ + 2C,C,Ck (4) 



or their connected counterparts 



and it is an easy exercise to find any particular quantity of tliis kind, for particular low values oi ii, . . . ,ik 

C2{N) = N^ 

Ci{N) ^2N^ + N, 

Ce{N) ^ hN^ + lOAT^, 

(5) 
or, say, for particular low values of TV: 

C'fci,...,*;„(l) = (fci + • • • + ^m — 1)!! if fci + . . . + fcm IS cvcu, and otherwise (6) 

The fact that each connected correlator is a non-trivial polynomial in TV is encoded in the idea of a genus expansion: 
different powers of N come from the fat-graph Feynman diagrams of different topology. The real problem is to 
find a closed expression for generic C and C"^°"" or, at least, for contributions of particular genera. 

As usual, straightforward approach is to construct a generating function and evaluate it exactly. Also as usual, 
resolvability of this problem depends on the clever choice of the generating function. The simplest choice is to 
consider resolvents 

p{zi, . . . , Zk) = = > -[-ri iT+T = (*•■ 7---t^ 7/ v) 

Z . ^^ ^ z^+ . ..z^^' \ zi-4> Zk-(p/ 



or 






p— (z„...,z,) = v(z,)...v(z,)iogz^ y: ^jv";;.;r - U^^it^---^^;;^)) («) 



These generating functions are very well adjusted to the use of Ward indentities for Gaussian matrix integrals, 
also known under the names of loop equations or Virasoro constraints, see [1] for a comprehensive review. The 
problem is, however, that only contributions of particular genera can be calculated in a closed form this way [2]: 

^{zi - Z2Y \yN[zi)yN[z2) 

^conn/, , , ^ _ 2JV(2:iZ2 + Z1Z3 + Z2Z3 + AN) 

Pa [zi,Z2,Z3) — — — — — [\i) 

yN(ziyyN(z2ryN[z3y 

where y%{z) — z^ — AN is what is called a spectral curve of the (in this case, Gaussian) matrix model. Given 
a spectral curve, Virasoro constraints can be re-interpreted as "topological recursion" rules [?, ?, ?], allowing to 
build pg(zi, . . . , Zn) one after another, recursively in the numbers of handles and punctures g and n. 

As one calculates more and more resolvents of particular genera, it starts getting clear that formulas become 
increasingly complicated: it is rather hard to imagine that generic expression for p or p^"^"^ can be obtained in any 
closed form. It turns out that a minor modification of the generating function makes it calculable! Moreover, these 
Harer-Zagier functions turn to be elementary functions. In [3] actually only a 1-point function was introduced 

and the difference from p(z) is the peculiar bifactorial in the denominator. As shown in [4], this result is by 
no means limited to a 1-point function: all Harer-Zagier functions 0(zi, . . . , z„) look comprehensible! To make 
expressions really elementary one should also get rid of N, by making additional generating function [2, 4]: 



oo 



1/1 1 \ A 1 



*WA, . i;^,WA.,A» = ^ y^-^ - —J ^ (^^71 . _,_(.,,),. (11) 



Then, for example [4], 

<^odd(.i,.2|A) = X.^A ^^l.^(2n + l)!!(2z. + l)!! ^-^+^^-^+^(^) " (1 - A)3/ ^ 1 - A + (1 + A)(z,^ + z|) ^''^ 

and so on, see s.2 for more examples. The purpose of this paper is to shed new light to calculability of such 
Harer-Zagier functions and make a step towards derivation of closed and generic expressions for all of them. The 
idea is to use a still another, third, type of the generating function, 

e(si,...,Sfc)= (tre^i^...tre^'=^) (13) 

and 

e^°""(si,...,Sfe)= ((tre^i^..tre^-'*)) (14) 

They turn out to be closely related to (Generalized) Kontsevich matrix model (GKM) [5, 1], which in Gaussian 
case possesses a remarkable duality [6] with the Gaussian matrix model, considered in the present paper. This 
duality led E.Breain and S.Hikami to a wonderful representation of e(s|A'^) in the contour integral form [7]: 

e(si,...,Sfc|iV)-TT— idM,e"'^' 1 + ^ det ■ ^ (15) 

fJl Si J \ UiJ l<l,J<k Ui - [Uj + Sj) 

The last determinantal factor here can be also rewritten as 

det ^ ^^{u,-u,){u.-u,+s,-s,) ^^g^ 

l<ij<fe Ui — (Uj + Sj) ^^. [Ui - Uj + Si){Ui — Uj — Sj) 

what makes the whole expression similar to LMNS integral [8] , which recently attracted a lot of new attention 
in the context of AGT-relations [9]: as already emphasized in [10] duality of [6, 7] should play a prominent 
role in the understanding of AGT relations, at least on the lines of the Dotsenko-Fateev-matrix-model approach 
[11]. Remarkably, eq.(15) possesses a direct generalization at least to cubic Kontsevich model [12], and questions 
now arise about further generalizations to arbitrary GKM and about the possible AGT-like applications of these 
generalizations. 

In this paper, however, we do not proceed in these very appealing directions. Instead, we concentrate on the 
application of (15) to evaluation of Harer-Zagier functions, which - in variance with p{z) and e(s) - are elementary 
fully calculable functions. Application of (15) for this purpose will be absolutely straightforward. First, conversion 
from A^ to A substitutes the factor 



n^ 



i=l 



N 
Si 



Ui 



by 



k 



1-An 1 



i=l 



Second, contour integrals just pick up contributions from poles of the simple rational function of Ui - in a 
very similar way to Nekrasov's celebrated calculation [13]. Third, the Harer-Zagier bifactorials are made from 
factorials in e**" by shifting Si — > Siti and introducing additional semi-integer-F-function factors with the help of 
i-integrations J^ e^*»/^(iii. 

In this way we reproduce the 1-point and 2-point Harer-Zagier functions from the Brezin-Hikami integrals. 
Generalization to 3-point and higher functions should be straightforward, but a clever method needs to be invented 
to convert the resulting multiple integrals into elementary functions (which, as we conjectured in [4], are always 
elementary functions of arctangent type). This conversion is trivial in the 1-point case; in the 2-point case for the 
purpose of such conversion we use scries representations, which generalize the elementary formula 



E79rrW7^-/„ /„ e-<''''sH2ct,t,)dt,dt, (17) 



^'^'^^^"T^fe? ^ 4fe(fc!)^ c^fe+ 1 

V^b^^ ^{2k + iy.{aby 

see s.3.2. It is not yet clear if this series approach is really promissing for 3-point and higher cases. Probably, some 
new ideas are required to reveal the general structure of fc-point Harer-Zagier elementary functions for fc > 3. 



2 Harer-Zagier formulas: a reminder 

Harer-Zagier functions (multidensities) in the finite N Gaussian model are defined as 

cl>{s^,...,su\N)^ Y, Cl°^X{N)^ -^ (f8) 

. ^ n{ii) n{ik) 

(we will consider only connected Harer-Zagier functions here) where the bifactorial normalisation is 

(■\ f (* ^ 1)-; ■* = even , , 

"W = | ^\\' ^= odd (19) 

As emphasized in the previous section, it turns out to be convenient to take a generating function with respect 
to the size of a matrix A'^ as well, converting it into a continuous parameter A: 

oo 

<^(si, . . . , Sk\\) = Y. ^'^<^(«i' • • • ' 'k\N) (20) 

Technically, this conversion is done with the following well-known finite and infinite summation formulas: since 

1 iN + k)\ ^ 

(f-A)fe+i ' ^ Nlkl ^ ^ ' 



N 



we have: 



N 





V ^'-^ 


N 


AT\N ^ 


(1-A)^ 


>> 

N 


2,N A(f + A) 
(1-A)3 


V^X^ 


A(1 + 4A + A2) 




(1-A)4 


A(1 + 11A+11A2 + A3) 






(22) 



Similarly, from an obvious formula 



E, , , (n + l)n(n — l)(n — 2) . . . (n — m) 
j(j - l)(j - 2) . . . (j - m) = ^ '-^ ^^— -^ ^ '- (23) 

(l.h.s. is a polynomial of degree to + f in n, which obviously vanishes at n = 0, 1, 2, . . . , m) it follows that 

n{n + 1) 



E^- 



2 



^ . n(n +l)(2n+l) 



6 

v^ .4 n(n + l)(2n+l)(3n2 + 3n-l) 

Ej =^ ^^ 35^ -■■■ (24) 



2.1 1-point function 

The 1-point Harcr-Zagier density has a form [2, 3, 4] 

This simple and beautiful formula admits generalizations to 2-point and higher functions, see [4] and below. 

2.2 2-point function 

The 2-point Harcr-Zagier density has a form 

4>{S1,S2\X) = 4>evenisi,S2\X) + 0odd(si , S2 | A) (26) 

where a convenient separation of the even and odd parts is introduced, and 

arctan ^i^=2VT^ 



(Podd[Si,S2X) = rr^T^ r- ^ ^_, ^^^ „ „^ 27 



siS2 { d d 

dsi ds2 



(i>even{si,S2\\) = — 3 ( ^^~„ '*2^— ) 4>odd{si, S2\\) (28) 



These formulas were derived in [4] with the help of integrable Toda equations. In the present paper we rederive 
them, using a different method (Brezin-Hikami contour integrals) see below. 

2.3 3-point function 

The 3-point Harcr-Zagier density has a form 

0(S1,S2,S3|A) = 0e-uen(si,S2,S3|A) + cj)mix{si ^ S2 , S^\X) (29) 

where (t>even{si, S2, ssl-^) is a totally even (even w.r.t each of its arguments) part of the function. We present the 
following expression for this totally even part (without a proof, at conjectural level): 



arctan ' ' 



/(r I si,S2,S3) 



1-Si 1 + S2+S3-Si\ I Y 1+t2.Si(.S2 + S3) 



(1+t2.Si(.S2 + S3))2 1+t2si(s2 + S3); / 52^3(1 -Si) 



1 +t2si(s2 + S3) 



( / SiS3(1-S2)t2 

arctan ' ' 



I-S2 I + S1+S3- S2\ I Y l+T2s2(si-hS3 



(1+t2s2(si+S3))2 l+r2s2(si+S3); / SiS3{l-S2) 



1 +t2s2(si + S3) 



/ / SiS2(l-S3)r2 

arctan ' ' 



I-S3 1 + Si+S2-S3\ I Y -^+^^'^3(S1+S2 



(1+t2s3(si+S2))2 l+r2s3(si +S2)/ / SiS2{l - S3)t^ 

Y 1+t2.S3(si+S2) 



Sl{s2 + 33) S2{si +53) 53(^1 +^2) 

S2S3(1 +t2.Si(s2 +S3)) SlS3(l +t2s2(si +S3)) Si S2 (1 + T2.S3 (si + S2)) 

(S1+S2)(S1+S3)(S2+S3) 



SlS2S3(l + T^(SIS2 + S1S3 + S2S3 - S1S2S3)) 

with a notation 

fir \ 81,32,33) = 5- — — (f>even {T3i,T32,Ts3 ———] , T = T\J\ - Si - S2 - S3 (31) 

' \-T^ S1S2S3 V 1+T/ 

This formula was not presented in [4], but it can also be deduced by the method of computer experiment, exploited 
in that paper, which becomes today a powerful approach to derivation of important results and conjectures in 
mathematical physics. It would be interesting to prove (derive) this formula, by using either Toda equations 
or Brezin-Hikami integrals. It is even more interesting to generalize to further generalize it to the 4-point and 
higher functions, and check the appealing conjecture that all the Harer-Zagier functions are expressed through 
actangents, as linear combinations, reminiscent of the Okounkov's cq. (??), which we remind in the Appendix. 

3 Harer-Zagier formulas from Brezin-Hikami integrals 

3.1 1-point function 

In the case of a 1-point function (15) is especially easy to derive [7], however already in this case the calculation 
is actually done in Gaussian Kontsevich model, not in Hermitian model itself - thus emphasizing the role of 
Hermitian-Kontsevich duality [6]. 

t / cA _1-vA2/9 fv AA-. , , X ^ / ^A-.. * 9/ ,N Clet 6 ^ ■' T— r —(h^j2ll 



Jn>cn ^J ''^' Ai<p)A{A) \^ 

^5:/A(0)e^^^^(-)-n^-^^/^^^'"-'^<^. . i-5:^e-^<. (32) 



where Vat = ycAsu{N) is the volume of unitary group and A is the matrix with eigenvalues 

if ^ - A, + 3Sk^ (33) 

This is an elementary essentially rational expression, however, it looks singular in the limit of yl — > 0. The limit 
can be easily taken if the formula is rewritten in a somewhat more complicated form of a contour integral [7] : 



y- TT a» - flj + s((5,fc - 5jk) T-j- g(a,+s5,fc)V2 ^ gsV2 TT ga?/2 V^ ^sa, TT aj - Clk - 3 

k i^k 

.sV2 -Q ga?/2 / ^us^^ -Q U-a., + 3 ^^^^ 

J . U CLi 



, .^. Qi — aj ^^ -^.-^ ^ -^J;- Qi — ak 

k i<j ■' I I k t^k 

= 6" 



For A = this turns into 



-e'"/^ Ie'''(^l + -Y du (35) 

S.N - 



Following [4] one can now pass to a generating function in N: f{N) -^ /(A) = ^^ X f{N). This could be 

Ni. 



done in different ways, for example, one could write ^^ J^fi^) instead. Our choice, however, makes calculation 



of contour integrals most convenient and elementary. Indeed, instead of (35) wc obtain: 



e(s|A) 



1-Mi + S) 



g 1-A 2 



(1-A) 



(36) 



inhere contour integral picks up a contribution from a single pole at u = t^ 



Harer-Zagier function differs by the substitution 



2k 



Since 



we have 



1 



2*fc! 1 



(2/c-l)!! {2k)l i2k)\Jo 



^2fe+lg-tV2rfi 



(s|A) - / e{.st\X)e-*"/Hdt = ^— ^ H e'^''('-^'')tdt = 



(i-xrJo 



A 


(1-A)(l- 


-A-(1 + A)s2) 



(37) 



(38) 



(39) 



what successfully reproduces (11). 

Note that in this case both exponential and Harer-Zagier generating functions are elementary. In fact from 
(36) one can also obtain an explicit formula for the 1-point Gaussian resolvent: 



p(z|A) = (tr 



z — (j) 



e(s|A)e -''ds = 

(1 - ^) Jo 



gl-A 2 ^^ds = 



iX 



(l-A)yr^^A2 \ Vl + A 



erf I iz\ll r I (40) 



We need here an expansion of error-function 



erf(z) 



e 2 ^^dx = 62/ e 2 dx 



x^=z^+t 1 r e-'dt 



^l 



2t 



(41) 



at large values of its argument, where it becomes a series in inverse powers of z. Adequate for working out this 
asymptotics is the last representation, which immediately implies: 



p{z\\) = 


^ A(l + A)'= (2A:-1)!! 


fe=0 ^ ^ 



(42) 



Using (22) one can easily check the consistency of this formula with the well known first terms of the p-expansion: 



P(^) 



< tr / > <tT(h^ > < tr (/)'' > 



N N^ 2N^ + N 



(43) 



It is also clear from (42) that the resolvent, including the contributions of all genera is represented by asymptotic, 
nowhere convergent series, and is not expressed through elementary functions. Expressible are particular terms 
p(p\^)(^z) of genus expansion [2], 



P(^)- 



z-y{z) N 2lN{z^ + N) lliV(135z4 + SSSiVz^ + ISSA^^) 



y'{z) y'Hz) 



y^^z) 



N N^ 2N^ + N , , 



with t/^(z) ~ z^ — AN. Particular terms of the genus expansion (44) can be directly deduced from the general 



formula (42). For this it should be re-expanded in monomials (1 — A) ^: 

.r-IAl - V ^(1 + ^)' (2fc - 1)!! _ ^ ^ (2fc - 1)!! 2^ f (1 - 0(1 - e/2)'- ,, ,, _ 

Py^\'^) Z^ n _ )^\k+2 ^2k+l 1^2-^ ^2k+l (l-\)s+lj ^fe+2 ^^ 

(21) .-^ 2fe(2fc - 1)!! {N + s)\ ^^ ni- g)(l - ^/2)\ ^ (21) ^ 2'^{2k-iy.l ^^ r (l-e/2)fc 

Z^ 22fe+l Z^Z^ M\g] '^ f ffe+2 ^ "^ Z^ ^2fe+l Z^^ f (-1 _ f\Ar^fc+2'^^ 

k s TV ■ ■ -^0 ^ fe ^ "'0 V s; S 



^.^^^^ 2fc(2fc-l)!! ,+i (l-e/2)^ 
Z^ i^z2fe+i(fc + l)! « (1-fW 



^ , - Kfc + 1)! ^ (1-ey 



«=o 



>r^ .^>r^ 2'=(2fc- 1)!! / iV'^+i iV'=-i iV'^-^j-g^ _ 2) N''-^{35P-77k+12) 

1^^ 1^- 



N 



^2k+i \[k + iy. 12(fc-2)! 1440(fc-4)! 128 • 81 • 35 • (fc - 6)! 

_^^^ 2'=^(2fc- 1)!! f N^+^ ^ iV*^-! ^ ^ N''-^ ^ iV'^-s 



^2fe+i [(fc + 1)! 12(fc-2)! \288{k-5y. 80(fc-4)! 

^fe-5 ]\[k-5 J\lk-5 



128 -81 -(A; -8)! 960 -(fc- 7)! 64 • 7 • (fc - 



In this manipulation eq.(21) is used twice: first in order to convert (1 — A)^*^^^ into a power series in A and second 
to find a sum over s, which is equal to (1 — £^)~ . The integration variable ^ is similar to 1 — A. 

Starting from 2-point functions, only Harer-Zagier functions are expressible in terms of elementary functions. 
Like more general Okounkov's formulas for cubic Kontsevich model, expressions for exponential correlators remain 
well-defined (convergent), but can be reduced only to integrals of elementary functions. Transition from resolvents 
to convergent exponential correlators is nothing but the application of the standard Pade summation method, 
while transition to Harer-Zagier correlators is a fancy modification of Pade method, leading - at least in the case 
of Gaussian correlators - not only to convergent, but to elementary expressions. In fact transition between the 
Harer-Zagier and exponential functions is also an ordinary Pade transform, only with a twice smaller parameter: 

/5(z) = ^^(2fc-l)!!B'= ^ ^(,)^A^S'= = ^ ^ e(,s) = A^ |^ ^ Ae^/^ (46) 

with A == 7]-Tx2), B = jY^s^ and z = s~^ 

3.2 2-point function 

Expression (47) remains almost the same for arbitrary exponential correlators, only for a correlator of n expo- 
nentials the label fc at the r.h.s. is a vector of the size n and 

e(.i,...,.„|7V) = -^ / ntre-^e--^V2,tr^,^^^ 1 ^^CA^^tri^,, ^^^^ 

and 

^ n 

if) =A + J2 So.5^k^ (48) 

Contour-integral representation, however, becomes a little more involved [7]: 

e(si,...,s„|7V)= I I rf) dwae"""" 1 + — M^ ^ --^ ^, (49) 

Q — 1 a<p r- ^ / 

SO that 

e(.i,...,.„|A) = n ^ /'^-.e--n TfTT ^ n i""""^^!""r^^'^"'^i (50) 



In the case of a two-point function, n = 2 the integral is two-fold: 



e(si,S2|A) 



= (s?+«2)/2 



S1S2 



(u — v){u — V + Si — 82) 



(1 — X)UV — \U32 — Af Si — \S1S2 



[u- V + si){u -V - 82) 



(51) 



and the main difficulty conies from quadratic factor in the denominator. It is actually linear both in u and v, 
thus the first contour integral, say, over v is simple: it picks up contributions from the three poles w = m + si, 
V = u — si and v = ii-X)u-\s • However, the last of the three gives rise to a new pole in m, and actually this 
is the only pole which contributes to the second integral over u: contributions of other poles in u cancel to zero. 
After all these intermediate calculations, the resulting expression for the exponential 2-point function is 



e(si,S2|A) = 



siA2 



[i-xy 



du /(l + A)(s2 + s2) Xs^sl 



«=o 



2(1 -A) u(l-A)2 



m(1-A)(si-S2) M^(l-A)' 



AsiS2 



1 _ m(1-A)(si-S2) 
AsiS2 



AsiS2 



and it is not so easy to evaluate, because u^^ is present in the exponent. The problem is slightly simplified by 
taking the odd part of the function, eodd{sl^ S2IA), which takes a somewhat more concise form 

eodd(si,S2|A) = -ie{si, S2\X) - e(si,-.S2|A) - e(-si,S2|A) + e(-si, -S2|A)j = 



e 2(i-A) 



du- 



Xs2{u'{l-Xf + Xsl) 



r„.o {X'^sl-u^{l-XY){sl-u^{l-XY) 
Rescaling then u 1— >■ u/ si and denoting S = siS2/(l — A), one finds 



sinh 



AS1S2 



u{i-xy 



X (i+A)(.;+4) 
eodd{si,S2\X) = - re ^d-^) 



1-A 






Using the simple identity 



we obtain a series 



S{u^ + XS^) 



(A2S'2-'u2)(5'2-'u2) ^u2'=+2 A-1 



E 



C2fc-|-1 \2fc-l-l _ 1 



SiU 



eodd(si,S2|A) 



X (i+A)(.? + e|) °f^ ~ ;^2fc+l _ I g2k+2m+l^^m 

e 2(i-A) \ \ 

'^ '^ A-1 m\(m + 2k + ly. 



1-A 



and this series can be used directly to calculate the Harer-Zagier function: 

00 00 
0odd(si,S2)= \ dtx dt2e^''^^^e^''^^^eodd{siti,S2t2\X} ^ 



00 00 



A ^ ^ \^M - 1 ^,„_^^^„. /^^\ ■■"■+'•"*' („, + 1)!(„. + ky 



T^Ei: 



k=Om=0 



A-1 



1-A 



k-\-m-\-l 



m\{m + 2k+iy. \l-i+^s2 l-i±fs2 



arctan 



yi^M(^f+i) 



(1-A)2 /l_l+A.2 



^l_l±A(,2+^,2) 



(52) 



In fact, the equality between the above double series and the 2-point arctangent function is not quite obvious and 
requires additional comments. A more natural series representation for the 2-point arctangent function is 



arctan 



^ab—c'^ 



\/ab — c2 



E 



{-fcy 



2k+l 



(2k + l){ab - c2)'=+i ^^ (06)*^+™+! 



-1 E (a 



^2fe2+2m+l fjyi^^l^y ry'ik+l 



i\k\ 2fc + l 



(53) 



with compact notations 7 = (1 — A)/(l + A), a = 1 — 75^, 6 — 1 — 7S2, c = S1S2. Comparing this natural double 
series with the above double series expansion of the contour integral, one observes a puzzling series identity 

^(p-fc)!b + fc + l)! ^ kl{p - k)l{2k + 1) plJo 

At A = 0, this identity becomes just a B-function evaluation 

while for generic A 7^ it can be considered as a hypergeometric identity and proven, say, by the systematic 
methods of Wilf and Zeilberger [14]. 

The even part of the contour integral can be treated similarly. As one can see, a quite complicated contour 
integral gives rise to a seemingly no-less-complicated series expansion, which - after the Harer-Zagier transform! 
- turns into an elementary function. However, we still do not posess any general expression for these elementary 
functions; even derivation of the simplest examples, as it is clear from the above calculation, is not quite easy. 
For higher-point functions, such calculations become increasingly more complicated, and remain a challenge. 

Appendix: Okounkov's formula for cubic Kontsevich model [12] 

Another important direction of generalization is from the Gaussian Kontsevich model to Generalized Kontsevich 
models of higher degree. In the case of cubic Kontsevich model 



f DX cxp ( -iTrX3 - J^TtAX^ o2fe+i 

7 r , T2k+i = — -— -TrA (56) 

JDX exp(-^Tr AX A 2fc + 1 



/i?Xexp(^-^_ 
where the conventional resolvents have a form 

°° 1 d 

PK{zi,...,Zm) ^ Mzi)---MZm)^OgZK _, A{z) :^ ^ »+3/2 ^ (^'^) 



n=0 



+3/2 dT2n+l 



a contour integral representation was obtained by Okounkov (see [12] for more details) for the exponential functions 
eii:(a;i, . . . ,Xk), which are connected with the above resolvents by Laplace transform: 

PK{zi,...,Zk) ^2 [[dsie ' 6^(51,..., Sfe) (58) 

*=i 
This contour integral representation has a form 

eK{s) =£(s), 

6^(51, S2) = £{si,S2) - £{si)£{s2), 

enisi, 82,83) = £{81,82, S3) - £(si)£(s2, S3) - £(s2)£(si, S3) - £(s3)£(si, 2) + £(si)£(s2)£(s3), 
and so on, where 



(59) 



n J_„3 



612"° f ,^ ( (ta-ta+lY {ta+ta+l)Sa 



a=l 



a Jta>n 



^(^i'---'^") = ^5^r^^ll^^ / dtgcxpl (60) 



These formulas can be considered as direct generalisations of Brezin-Hikami contour integrals for the Gaussian 
Kontsevich model. It is puzzling to ask, whether similar integral representations exist for arbitrary monomial 
potentials in GKM, or maybe even for generic non-monomial potentials. Of interest is also to find analogues of 
elementary multidensities in these models. 
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